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We examine the consequences of the effective field theory (EFT) of dark matter-nucleon scattering 
for current and proposed direct detection experiments. Exclusion limits on EFT coupling constants 
computed using the optimum interval method are presented for SuperCDMS Soudan, CDMS II, and 
LUX, and the necessity of combining results from multiple experiments in order to determine dark 
matter parameters is discussed. We demonstrate that spectral differences between the standard 
dark matter model and a general EFT interaction can produce a bias when calculating exclusion 
limits and when developing signal models for likelihood and machine learning techniques. We also 
discuss the implications of the EFT for the next-generation (G2) direct detection experiments and 
point out regions of complementarity in the EFT parameter space. 


I. INTRODUCTION 


Astrophysical and cosmological evidence indicates that 
the majority of the matter in the universe takes the form 
of non-luminous particles called dark matter, though the 
exact nature of the dark matter particle remains un- 
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known [1]. A generic weakly-interacting massive parti¬ 
cle (WIMP) is a very attractive dark matter candidate 
[2]. Numerous experiments are engaged in efforts to de¬ 
tect rare collisions between WIMPs and target nuclei in 
terrestrial detectors. Results from DAM A [3], CoGeNT 
[4], CRESST-II [5], and CDMSII Si [6] can be inter¬ 
preted in terms of interactions of WIMPs with masses 
of 6-30 GeV/c^. A similar range of masses could also ac¬ 
count for a possible excess in the gamma-ray flux near the 
galactic center in Fermi-LAT data [7, 8]. Under standard 
assumptions for spin-independent WIMP-nucleon inter¬ 
actions, however, such interpretations are difficult to rec¬ 
oncile with the limits set by GDMSlite [9], SuperGDMS 
[10], LUX [11], and PIGO [12], 

Standard WIMP scattering calculations make simpli¬ 
fying assumptions about the type of interaction between 
the nucleon and the dark matter particle: typically only 
isospin-conserving spin-independent couplings, or spin- 
dependent couplings to either the proton or neutron are 
considered. This results in constraints on the three cor¬ 
responding WIMP-nucleon cross sections. Relaxing such 
assumptions can suppress the interaction for some tar¬ 
get elements by orders of magnitude relative to others 
[13]. In particular, assuming different spin-independent 
dark matter couplings to protons, /p, and neutrons, /„, 
can reconcile much of the tension between the GDMSII 
Si allowed region and the SuperGDMS Soudan and LUX 
exclusion limits [14]. However, such solutions often re¬ 
quire a high degree of fine-tuning. 

In addition, the calculation of dark matter scattering 
rates typically assumes a Maxwellian velocity distribu¬ 
tion [15]. As shown in [16, 17], N-body simulations are 
not well described by such a distribution. Gonsequently, 
alternate halo models have been proposed. One such ve¬ 
locity distribution is discussed in [18, 19] and takes the 
form 


erties in addition to the standard spin-independent and 
spin-dependent cases. It also explicitly includes isospin 
interference and interference between operators, creating 
a rich parameter space of possible dark matter interac¬ 
tions that are very sensitive to the specific choice of de¬ 
tector material. 

The EFT framework parametrizes the WIMP-nucleus 
interaction in terms of fourteen operators, Oi, which are 
listed in Eq. 2 and include the standard spin-independent 
and spin-dependent interactions. These operators feature 
explicit dependence on (the relative velocity between 
the incoming WIMP and the nucleon) and the momen¬ 
tum transfer (f, in addition to the WIMP and nucleon 
spins, and Sn- Note that O 2 is not considered since it 
cannot arise from the non-relativistic limit of a relativis¬ 
tic operator at leading order. In addition, each operator 
can independently couple to protons or neutrons. We for¬ 
mulate this isospin dependence in terms of isoscalar and 
isovector interactions, following the conventions of [22]. 
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for dark matter velocities smaller than the galactic escape 
velocity Vesc- For values of Vq/V(,sc and p consistent with 
N-body simulations, this function falls off faster than the 
standard Maxwellian distribution. This difference can 
significantly affect the expected dark matter event rate, 
especially for low-mass WIMPs for which experiments 
are only sensitive to the high-velocity tail of the distri¬ 
bution. It has been shown that choosing certain values 
for the parameters of this alternate halo model can rec¬ 
oncile the tension between GDMS II Si and XENONIOO 
[20], though it cannot also account for the tension with 
LUX because of that experiment’s lower energy thresh¬ 
old. 

Recently, an effective field theory (EFT) approach for 
WIMP scattering has been developed that considers all 
leading-order and next-to-leading order operators that 
can occur in the effective Lagrangian that describes the 
WIMP-nucleus interaction [21-23]. This formalism intro¬ 
duces new operators that rely on a range of nuclear prop¬ 
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These operators contribute to six types of nuclear re¬ 
sponse functions. The spin-independent response is de¬ 
noted by M and is typically the strongest of the six func¬ 
tions since it is related to the number of nucleons in the 
target nucleus. The main contribution to this response 
comes from the standard spin-independent operator C>i, 
but it also contains higher-order contributions from op¬ 
erators O 5 , Os, and On. There are two spin-dependent 
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responses, S' and S", that correspond to projections of 
spin parallel and perpendicular to the momentum trans¬ 
fer. A linear combination of these two responses yields 
O 4 , which is related to the standard spin-dependent re¬ 
sponse. Many of the other operators also appear in one of 
these two responses. A novel type of response introduced 
in the EFT, A, is related to the net angular momentum 
of an unpaired nucleon and contains contributions from 
operators O 5 and Os- A second novel response is 
which is sensitive to the product of angular momentum 
and spin. This response tends to favor heavier elements, 
and the most dominant contribution to this response is 
from O 3 . The last response considered in the EFT, 
contains contributions from operators O 3 , O 12 , and O 15 . 

is discussed less frequently in the literature since it is 
difficult to find a model that produces this response, but 
we consider it here for completeness. 

The EFT also includes two operator-operator interfer¬ 
ence terms: E'A and Md*". E' interferes with A because 
velocity-dependent responses are sensitive to properties 
such as angular momentum that depend on the motion 
of the nucleon within the nucleus. This interference term 
is particularly significant for germanium, which has large 
responses to both E' and A. The E'A response contains 
interference between O 4 and O 3 , as well as between Og 
and Og. In addition, since both M and $" are scalar 
responses, interference between the two can be signifi¬ 
cant, especially for elements like xenon that have large 
responses to both. The response contains interfer¬ 
ence between operators Ox and O 3 , operators On and 
O 12 , and operators On and O 15 . 

Since the various responses are related to different nu¬ 
clear properties, the strength of the resulting interaction 
can vary by many orders of magnitude. The expecta¬ 
tion values of these properties are listed in [21]. For in¬ 
stance, the spin-dependent responses E' and E" depend 
on the square of the spin of an unpaired nucleon, which 
ranges from 5 x 10“® for protons in germanium (which 
has one isotope with an unpaired nucleon, which is a 
neutron) to 0.2 for protons in fluorine (which has an un¬ 
paired proton). The angular momentum of a nucleon, 
which governs the strength of the A response, ranges 
from 0(1 X 10“^) to 0(1), while {L ■ 5)^, which gov¬ 
erns the strength of the $" response, ranges from 0.1 for 
light nuclei to several hundred for heavier nuclei. The 
strongest response is M, which is related to the square 
of the number of nucleons. 

The strength of an EFT interaction is parametrized 
by numerical coefficients, cj", associated with each op¬ 
erator Oi, where r = 0 or 1 denotes the isoscalar 
(c? = V 2 (cf + ^'i)) isovector {c] = ^/ 2 {c^ — c")) 
combinations, respectively. The coefficients have dimen¬ 
sions of 1 /energy^, so we multiply by the weak mass scale 
i'rriweak = 246.2 GeV) to produce dimensionless quanti¬ 
ties. The c[ are related by a change of basis to generalized 
versions of /„ and fp and can take on any value, posi¬ 
tive or negative. The coefficients appear as c[cj in the 
interaction, indicating that operators interfere pair-wise. 


at most. 

This paper discusses the EFT approach in the context 
of current and proposed direct detection experiments. 
We present exclusion limits on EFT operator coefficients 
using the optimum interval method [24, 25]. We discuss 
the differences in energy spectra that arise for arbitrary 
EFT interactions and examine how this energy depen¬ 
dence may affect future experiments if WIMP candidate 
events are observed. We also consider the variation in 
interaction strength across the elements commonly used 
as direct detection targets and discuss possible ways of 
exploring interference using experimental results. Fi¬ 
nally, we discuss the implications of this effective field 
theory for the next-generation (G2) direct detection ex¬ 
periments, SuperGDMS SNOLAB and LZ. 


II. EXCLUSION LIMITS ON A SET OF EFT 
OPERATORS 

The strength of the interaction in the EFT frame¬ 
work is governed by a set of 28 numerical coefficients 
corresponding to the 14 operators, one for each isospin. 
Others have attempted to find global fits in this multi¬ 
dimensional parameter space, combining data from many 
direct detection experiments [28]. Since the parameter 
space is large and relatively unconstrained by current 
experiments, we choose instead to calculate exclusion 
limits on the coefficients for individual EFT operators 
for three different target elements: germanium (Super- 
CDMS Soudan and CDMSII), silicon (CDMSII), and 
xenon (LUX). This paper presents the first EFT experi¬ 
mental result that includes all three target elements that 
will be used in the G2 experiments. 

We use the optimum interval method to calculate 90% 
upper confidence limits on the numerical coefficients of 
EFT operators. The optimum interval method incor¬ 
porates information about the candidate event energies 
and energy-dependent detection efficiencies, which can 
yield stronger exclusion limits in the presence of unknown 
backgrounds than likelihood methods that consider only 
a single energy bin in the presence of backgrounds. This 
is particularly important here because of the spectral dif¬ 
ferences that can arise from different EFT interactions. 
We consider a single operator at a time and present 
the exclusion limit on the square of the EFT coefficient, 
which is proportional to the total interaction cross sec¬ 
tion. We compare the effects of two halo models on the 
limits. The first uses standard halo assumptions as in 
[15], with a WIMP mass density po = 0.3 GeV/c^/cm^, 
most probable WIMP velocity of 220km/s, mean circular 
velocity of the Earth with respect to the galactic center 
of 232km/s, galactic escape velocity of 544km/s, and a 
velocity distribution that correctly takes into account the 
effect of the Earth’s velocity on the escape-velocity cut¬ 
off [29]. The second halo model uses the functional form 
of Eq. 1 with p = 2.7 and vg/vesc = 0.6875, determined 
by fits to the Eris simulation of a Milky-Way-like galaxy 
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FIG. 1. Upper limits on the dimensionless isoscalar coefficients C 3 (left) and Cs (right) as a function of WIMP mass for 
SuperCDMS Soudan (light blue) [10], CDMSII Ge reanalysis (dark blue) [? ], and CDMSII Si (red) [27], and estimated limits 
for LUX (black) [11], for the Maxwellian halo (solid) and an alternate halo model (dashed). 


[30], and other halo parameters as above. 

Figure 1 shows the upper limits for two example oper¬ 
ators, isoscalar operators O 3 (left) and Og (right), as a 
function of WIMP mass. Limits on all operators for a 
small range of masses can be found in Table I. Limits on 
all operators for a small range of masses can be found in 
Table 1. Solid lines correspond to the Maxwellian halo, 
whereas dashed lines show the limit calculated assum¬ 
ing the alternate velocity distribution function discussed 
above. The SuperCDMS Soudan, CDMSII Ge (reanal¬ 
ysis), and CDMSII Si limits use the candidate events, 
thresholds, and detection efficiencies discussed in [ 10 ], [? 
], and [27] respectively, while the estimated LUX limit 
assumes zero observed events and functional form for the 
detection efficiency that follows a hyperbolic tangent ver¬ 
sus energy centered at 2.5 keVnr but with a step function 
cutoff that goes to zero below SkeVnr- 

Because of the different nuclear responses for the three 
target elements considered, the relative strength of the 
limits varies from operator to operator. In particular. 
Os (Fig. 1, right) includes contributions from the A re¬ 
sponse, which is greater in germanium than in silicon or 
xenon. This contribution strengthens the SuperCDMS 
Soudan constraint relative to LUX and CDMSII Si. In 
addition, the shape of the curve for a single target ele¬ 
ment changes from operator to operator. For example. 
Os depends on the square of the momentum transfer, 
naturally suppressing the event rate at low energies. As 
a result, the limits at low WIMP mass for O 3 are weaker 
than for other operators. 

The difference between the two WIMP velocity distri¬ 
butions becomes apparent when the only events expected 
above the detection thresholds are due to WIMPs in the 
high-velocity tails. Since both CDMS and LUX have 
thresholds of a few keV, this disparity appears only at 
the lowest WIMP masses. The difference is also more 


pronounced for LUX, since its target nucleus, xenon, is 
heavier than silicon or germanium. A dark matter parti¬ 
cle must have a higher velocity to deposit a given recoil 
energy in xenon than in germanium or silicon; higher- 
energy recoils become comparatively rarer. For the Su¬ 
perCDMS Soudan result, the difference in velocity dis¬ 
tributions leads to a factor of two difference in the limit 
around 4GeV/c^, whereas for LUX, the difference in 
velocity distribution leads to a factor of two difference 
around 7 GeV/c^. 



_0 ,„2 
^1 ^weak 


FIG. 2. Polar limits on Oi isospin for SuperGDMS Soudan 
(blue) [10], LUX [11] (black), and CDMS II Si (red) [27] at a 
WIMP mass of 6 GeV/c^ 

Since the EFT explicitly includes isospin dependence, 
we can also use the optimum interval method to set polar 
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TABLE I. SuperCDMS and CDMSII 90% confidence level upper limits on the square of the dimensionless EFT coefficient for 
pure isoscalar interaction for a 10 GeV/c^ (300 GeV/c^) WIMP for all isoscalar EFT operators. The upper limits vary in 
accordance with the relative strength of the interaction in silicon and germanium. 


limits on isospin. For a given WIMP mass and a given 
angle between the isoscalar and isovector components of 
an operator, we set a 90% upper confidence limit on the 
isoscalar-isovector radius. Varying the polar angle pro¬ 
duces exclusion ellipses in the isoscalar-isovector plane, 
as in Fig. 2, which shows limits for operator Oi and a 
6 GeV/c^ WIMP. The major axis of each ellipse corre¬ 
sponds to the value of c^/cj that yields maximum sup¬ 
pression of the scattering rate. Note that although the 
exposures for CDMS and SuperCDMS are significantly 
lower than for LUX, there are regions of parameter space 
allowed by LUX but excluded by SuperCDMS and CDMS 
at 90% confidence. This example demonstrates that a 
combination of experiments using several target nuclei 
can constrain the EFT parameter space better than any 
single experiment. 

III. EFFECT OF EFT ENERGY DEPENDENCE 
ON STANDARD LIMITS 

Because of the additional momentum dependence of 
several of the EFT operators, the differential event rate 
for an arbitrary dark matter interaction could be very dif¬ 
ferent than for the standard calculation. Consequently, 
it is possible that a limit-setting algorithm that expects 
the (approximately) exponential event rate of the stan¬ 
dard spin-independent interaction could misinterpret a 
potential signal from a more general EFT interaction as 
background. 

To demonstrate the possible bias that could arise from 
assuming the standard spin-independent event rate when 
setting limits, we perform simulated experiments assum¬ 
ing that the dark matter scattering is purely due to a 
single isoscalar EFT operator. Figure 3 shows the co¬ 


added results of 100 simulated experiments sampled from 
the energy spectrum of isoscalar O 3 scattering in germa¬ 
nium for two different dark matter masses, assuming an 
energy-independent (or “flat”) detection efficiency. The 
operator coefficients were set to give each simulated ex¬ 
periment an expectation value of 10 events. This expec¬ 
tation was then convolved with a Poisson distribution to 
select the number of events for a given simulated exper¬ 
iment. 

Unlike the standard spin-independent event rate, the 
event rate for O 3 depends on the square of the momen¬ 
tum transfer, so the event rate is suppressed at low recoil 
energies. This effect is illustrated in Fig. 3, where the 
black, cyan, and magenta curves show the standard spin- 
independent scattering rate for a range of WIMP masses 
and the blue histogram corresponds to the simulated 
spectrum expected from O 3 interactions. For the case 
of a 10 CeV/c^ WIMP mass, the distribution of events 
is more closely matched by the spin-independent rate for 
a higher-mass WIMP. For the 300 CeV/case, no spin- 
independent rate calculation matches the observed spec¬ 
trum of events; if experimenters only consider the spin- 
independent WIMP rate, they may erroneously conclude 
that they have observed an unexpected background or 
incorrectly measured their detection efficiency as a func¬ 
tion of energy. 

We calculate the 90% confidence level upper limit on 
the spin-independent cross section for each simulated ex¬ 
periment using the optimum interval method [24, 25] and 
the standard Maxwellian halo model with halo parame¬ 
ters as above with no background subtraction. Each sim¬ 
ulated experiment was assumed to have an exposure of 
1000 kg days and a flat efficiency of 60% between 1 and 
lOOkeVnr- The distribution of limits is shown in Eig. 4. 
Eigure 4 also shows the resulting median limit from sim- 
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EFT event sample for 10 GeV WIMP In Ge 



EFT event sample for 300 GeV WIMP In Ge 
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FIG. 3. Co-added energy spectrum from 100 simulated experiments (blue histogram) assuming the dark matter interaction 
proceeds according to the isoscalar O 3 operator for a 10 GeV/(left) and a 300 GeV/WIMP (right). The detection efficiency 
is assumed to be independent of energy. The smooth cyan, magenta, and black curves show the expected spectrum for the 
standard spin-independent rate for several WIMP masses, while the dashed dark blue curve shows the CI 3 spectrum from which 
the simulated experiments were sampled. 



FIG. 4. Distribution of 90% confidence level upper limits calculated using the optimum interval method for the simulated 
experiments discussed in Sec. 3 and shown in Fig. 3, sampled from the event rate for isoscalar O 3 . Shaded blue bands show 
the 68 % and 95% confidence level uncertainty on the distribution. The zero-background Poisson limit is shown in magenta. 


ulated experiments sampled from the spin-independent 
distribution in black. 

The distribution of limits on the spin-independent 
cross section for the simulated experiments sampled 
from the O 3 energy spectrum deviates from the zero- 
background limit shown in magenta as well as from 
the mean limit derived from similar simulated experi¬ 
ments sampling from the spin-independent rate. As ex¬ 
pected, the simulated-experiment limits are weaker than 
the zero-background limits due to the presence of can¬ 
didate events. However, because the energy distribu¬ 
tion of the candidate events sampled from O 3 is differ¬ 
ent than the expected spin-independent rate, the limits 


also deviate from the expected shape for the true spin- 
independent experiment. 

In the 10 GeV/c^ case, we expect the limit to be weak¬ 
est around a mass of 10 GeV/c^, where the rate expected 
by the limit algorithm matches the observed event rate. 
However, because the observed events due to O 3 scatter¬ 
ing are skewed towards higher recoil energies, the limit 
tends to be weaker at larger WIMP masses where the 
tail of the spin-independent event rate extends to higher 
recoil energies. For the 300 GeV/case, the distribu¬ 
tion of limits agrees with the Poisson zero-background 
limit at low masses; the observed events occur at recoil 
energies that cannot be produced by a low-mass WIMP. 
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At higher masses, the distribution of limits is still close 
to the zero-background limit because the shape of the 
observed spectrum is very different from the expected 
spin-independent WIMP rate. 

The difference in the limits between the spin- 
independent and EFT cases demonstrates the impor¬ 
tance of correctly modeling the expected WIMP signal. 
Algorithms that assume the standard spin-independent 
rate when calculating limits will interpret events from 
EFT interactions with different spectral shapes as back¬ 
ground, and thus, this assumption could lead to a bias in 
the exclusion limits reported by experiments, especially 
in the case where events are observed. 


IV. INTERFERENCE IN THE EFT 
PARAMETER SPACE 

A. General interference framework 

The EFT framework also provides a more general de¬ 
scription of interference among operators such as the 
“xenophobic” isospin violation case discussed in the lit¬ 
erature [13]. It not only allows for interference between 
the isospin components of individual operators, but also 
among different operators. The generalized interference 
can be written as a matrix equation in the large EFT pa¬ 
rameter space, but because operators interfere in pairs, 
and only certain pairs interfere, this large matrix can be 
decomposed into block-diagonal form. We consider the 
2 x 2 case of isospin interference and the 4x4 case of 
isospin and operator-operator interference. 

The generalized amplitude for the 4x4 case can be 
written as the product of the vector of operator coeffi¬ 
cients c[ with the amplitude matrix, where superscript 0 
and 1 indicate isoscalar and isovector, respectively, and 
the subscripts indicate the operator being considered: 


400 401 400 401 

^ij ^ij 


'c° ' 

410 411 410 411 

^ij ^ij 


G 

400 401 400 401 

ji ji jj jj 


c° 

^3 

410 411 410 411 

L ^3j ^33 J 


cl 

L ^ J 


The amplitudes AIJ are the product of the WIMP and 
nuclear response functions for the interaction specified 
by c[ and cj and depend on properties such as tar¬ 
get element, WIMP mass, WIMP spin, WIMP velocity, 
and nuclear recoil energy. We evaluate the AIJ with¬ 
out integrating over the dark matter velocity distribu¬ 
tion to avoid introducing more variables. Amplitudes are 
summed over the isotopes for a given element according 
to their natural abundances. 

Finding the eigenvectors of this matrix will give the 
“principal components” of the interaction space. We ex¬ 
pect that three of the four eigenvalues should be small, 
since the matrix for a single isotope is an outer product 


and therefore should have a single nonzero eigenvalue. 
The vector with the largest eigenvalue corresponds to 
the maximal amplitude for scattering in the interference 
space under consideration, while the three small eigen¬ 
values correspond to local extrema in the scattering am¬ 
plitude which tend to suppress the event rate. To be 
maximally sensitive to the parameter space for a given 
interference case, we would like to choose target ele¬ 
ments whose constructive interference eigenvectors span 
the space of interactions. 

As an example, we first consider isospin interference for 
a single operator in an already well-understood case. Fig¬ 
ure 5 shows the constructive isospin interference eigen¬ 
vectors for scattering via operator O 4 (the standard spin- 
dependent operator) for several elemental targets, evalu¬ 
ated at a WIMP mass of lOOGeV/c^ and nuclear recoil 
energy of 100 keV. The vectors are plotted in the space 
of the isoscalar coefficient versus the isovector coefficient. 
The proton-neutron space can be recovered from this ba¬ 
sis via a 45-degree rotation. The amplitude in a given 
direction indicates the target’s response to that operator 
and illustrates the sensitivity of each material to the cor¬ 
responding operator. In addition, if we were to plot polar 
limits as in Fig. 2 for O 4 , we would see that the direc¬ 
tion of the constructive interference vector corresponds 
to the minor axis of the ellipse. In the two-dimensional 
case, the destructive interference vector is perpendicular 
to the constructive vector and corresponds to the major 
axis of the ellipse in a polar limit plot. 



FIG. 5. Constructive interference eigenvectors for 2D O4, 
isospin interference. Proton-dominated interactions occur 
along the x = y diagonal, while neutron-dominated inter¬ 
actions occur along the x = —y diagonal. 


Since O 4 is the standard spin-dependent operator, we 













see that the constructive interference eigenvectors fall 
into two categories based on the nucleon content of the 
target nucleus. The elements with unpaired protons (flu¬ 
orine, sodium, and iodine) have maximal scattering rates 
when the interaction is proton-dominated, correspond¬ 
ing to c° = c^. On the other hand, the elements with 
unpaired neutrons (germinum, xenon, and silicon) have 
maximal scattering rates when the interaction is neutron- 
dominated, corresponding to c° = —c^. Consequently, 
to span this space and therefore be maximally sensi¬ 
tive to all possible spin-dependent interactions, we should 
choose one element each from the neutron- and proton- 
dominated sets. 

We can apply this same procedure to the more general 
4D case to demonstrate the complementarity of the dif¬ 
ferent target elements. As an example. Figure 6 shows 
all 2D projections of the four-dimensional eigenvectors 
in the interference space for Os and Og, evaluated for a 
WIMP mass of lOOGeV/c^ and nuclear recoil energy of 
30 keV. The eigenvectors for scattering in silicon, germa¬ 
nium, xenon, iodine, and sodium indicate that they are 
most sensitive to various combinations of isoscalar and 
isovector Os scattering. However, the vector for fluorine 
shows that it is sensitive to both Os and Og. This varia¬ 
tion across targets allows different experiments to probe 
different regions of the EFT parameter space, increasing 
the overall sensitivity of the direct detection method. 



FIG. 6. Constructive interference eigenvectors for 4D OgjOg 
interference. 


To demonstrate the effect of this four-dimensional in¬ 
terference on the differential event rate, we evaluate the 
event rate using the operator coefficients from two four¬ 
dimensional interference eigenvectors from Fig. 6 that 


point in different directions in the parameter space. Fig¬ 
ure 7 shows the differential event rate for several targets 
evaluated at the constructive interference vectors for flu¬ 
orine (top) and germanium (bottom) for Os/Og interfer¬ 
ence. Since the fluorine eigenvector is not parallel to the 
germanium eigenvector, the germanium event rate eval¬ 
uated at the fluorine vector is suppressed and vice versa. 
In addition, since the xenon and germanium eigenvectors 
are nearly parallel in this case, the two event rates are 
comparable at the 30 keV nuclear recoil energy at which 
the eigenvectors are evaluated. 




FIG. 7. Differential event rate evaluated at the OsjOg con¬ 
structive interference vector from Fig. 6 for fluorine (top) and 
germanium (bottom). 

This example shows the large variation in signal 
strength that can occur for different combinations of op¬ 
erators. In this case, varying the coefficients from the 
germanium eigenvector to the fluorine eigenvector led to 
an order-of-magnitude suppression of the rate in germa¬ 
nium, silicon, and xenon, and a change in the energy 
spectrum for fluorine. Similar suppression can also occur 
for the other interference terms present in the effective 
field theory. 
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B. Comparison of target elements for G2 direct 
detection experiments 

Three target elements will be used in the upcoming G2 
experiments: germanium, silicon, and xenon. Under the 
standard spin-independent scattering framework, where 
the rate scales as experiments that use xenon as 

a target element have the greatest sensitivity for WIMP 
masses above a few GeV/c^. However, in order to probe 
operators dependent on other nuclear properties, the 
complementarity of the three G2 target elements merits 
further investigation. 

When considering the possible observations the G2 
experiments may make, the difference in experimental 
parameters such as detector mass and trigger threshold 
must also be taken into account. The proposed LZ de¬ 
tector will have a 5600 kg fiducial mass of xenon, while 
SuperGDMS will be operating 57 kg of germanium and 
silicon. Figure 8 shows the relative event rates for the 
three G2 target elements assuming scattering proceeds 
via a single isoscalar EFT operator. This figure only 
shows the relative WIMP rates for the G2 experiments; 
background rates are not taken into consideration. Note 
that the true interaction, which may come from a lin¬ 
ear combination of operators, could enhance or suppress 
these rates. 

We normalize the event rate so that SuperGDMS 
Ge observes one event for a given operator in the 
3GeV/c^ case and LZ observes one event for a given 
operator in the 10 and 300 GeV/cases. The LZ rate 
(black) assumes a 5600 kg hducial mass, an exposure of 
1000 days, a 100% trigger efficiency between 1 and 30 
keVnr, and a flat 50% nuclear-recoil selection efficiency. 
The SuperGDMS Ge rate (blue) assumes 50 kg of germa¬ 
nium operating in standard iZIP mode [31], an exposure 
of 1000 days, and a 100% trigger efficiency between 0.5 
and lOOkeVnr, and a flat 60% combined fiducial-volume 
and nuclear-recoil selection efficiency. The SuperGDMS 
Si rate (red) assumes 1 kg of silicon and an exposure of 
1000 days. Since the silicon detectors will be operated in 
high-voltage mode [9] , the trigger threshold will be much 
lower, so we assume a 60% combined trigger and fiducial- 
volume efficiency up to 50 keVnr, with a trigger threshold 
of 70 eV. We also plot the event rate for SuperGDMS Ge 
high-voltage (light blue) for the 3GeV/c^ WIMP case. 
For the SuperGDMS Ge high-voltage detectors, we as¬ 
sume a target mass of 6 kg, trigger threshold of 80 eV, and 
all other parameters identical to SuperGDMS Si high- 
voltage. 

Though silicon, germanium, and xenon have similar 
nuclear properties (e.g., all three have isotopes with 
unpaired neutrons), the variation in the event rate 
across operators and target elements is large. For the 
3GeV/c^ case, the strength of the silicon signal rela¬ 
tive to the germanium signal varies by three orders of 
magnitude, depending on the operator assumed. The 
signal in LZ is very close to zero for such a low-mass 
dark matter particle because the velocity required for a 


3GeV/c^ WIMP to deposit energy above the assumed 
1 keV threshold is greater than the galactic escape ve¬ 
locity. However, for WIMP masses above a few GeV/c^, 
LZ’s exposure, which is approximately 100 times larger, 
leads to event rates that are enhanced by approximately 
the same factor. In addition, the relative rate for Super¬ 
GDMS Si HV becomes smaller at higher masses, since, 
by design, it is mainly sensitive to the small energy de¬ 
positions produced by low-mass WIMPs. 

To examine the effects of the different possible inter¬ 
actions for experiments with similar fiducial masses, we 
also plot the event rate per time per target mass (Fig. 8 , 
bottom). Here, we see that both Ge and Si SuperGDMS 
detectors operating in high-voltage mode are more sen¬ 
sitive to low-mass WIMPs because of their lower thresh¬ 
olds. In particular, the germanium high-voltage rate per 
kg day (light blue) is nearly an order of magnitude larger 
than the standard germanium iZIP rate (blue) for certain 
operators. For higher masses, the rates for xenon (black) 
and germanium are comparable within an order of mag¬ 
nitude, but the nuclear properties of silicon (red) make 
it less sensitive to these interactions. In addition, Su¬ 
perGDMS Ge sees a modest enhancement to the overall 
event rate at high WIMP masses where the distribution 
of events extends beyond the assumed 30 keVnr upper 
limit for LZ. This effect is most prominent for opera¬ 
tors such as O 3 and O 15 , which have a dependence 
that suppresses the rate at low energies, though it is not 
enough to overcome the effects of LZ’s larger target mass 
in the total number of events. 

The variation in signal strength across target elements 
in this effective field theory solidihes the case for us¬ 
ing multiple targets to detect dark matter. Previous 
work has shown that complementary target elements 
can break the degeneracy between the standard spin- 
dependent and spin-independent operators [32, 33], and 
others have shown that this concept can also be applied 
to the larger EFT parameter space [28]. Such consid¬ 
erations are particularly important when incorporating 
the effect of interference on the event rate. Because of 
the presence of both isospin interference and operator- 
operator interference, there are many combinations of 
interactions that may greatly suppress the event rate for 
one particular element. Even if a single experiment sees 
no signal due to interference effects, a complementary 
target with different nuclear properties may still observe 
events. 

To demonstrate the effect of interference on the rela¬ 
tive event rate, we determine regions of extremal interfer¬ 
ence in germanium using the principal component anal¬ 
ysis method detailed above. The event rate suppression 
relative to Oi for the three G2 experiments for germa¬ 
nium constructive interference and destructive interfer¬ 
ence are shown in Fig. 9 for WIMPs with masses of 3, 
10, or 300 GeV/c^, assuming the standard Maxwellian 
halo model and the same experimental parameters as in 
Fig. 8 . Again, this figure does not consider the relative 
background rates for the three experiments. We consider 
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SGeVWIMP lOGeVWIMP 300 GeV WIMP 



FIG. 8. Relative event rates for LZ (black), SuperCDMS SNOLAB Ge iZIP (blue), and SuperCDMS SNOLAB Si (red), 
normalized to 1 observed event in SuperGDMS Ge (3GeV/c^) or LZ (10, 300GeV/c^). From left to right are shown the 
rates for a 3, 10, and 300 GeV/WIMP, assuming isoscalar interactions and the standard Maxwellian halo model. The 
3GeV/c^ case also shows the rates from SuperGDMS SNOLAB Ge high-voltage (light blue), which has similar parameters to 
SuperGDMS Si high-voltage, but a target mass of 6 kg. The top row shows cumulative event rates, while the bottom row shows 
events per time per target mass. True interaction strengths may differ from this calculation since the interaction may proceed 
via a linear combination of operators. 


all seven possible cases of four-dimensional operator- 
operator interference. The sum of the squares of the EFT 
coefficients is equal for all cases presented; however, be¬ 
cause of the relative strength of various operators and 
the presence of interference, the rate can be suppressed 
by many orders of magnitude. We characterize the inter¬ 
ference using the magnitude of the eigenvalue: the largest 
eigenvalue corresponds to the maximally-enhanced event 
rate, while small eigenvalues correspond to varying levels 
of destructive interference. 

The relative event rates in Fig. 9 indicate that con¬ 
structive interference can only modestly enhance the 
event rate. In the case of O 1 /O 3 interference, the maxi¬ 
mal rate is only ^ 1.5% larger than the pure Oi rate. For 
operators such as O 4 that depend on the spin of a nucleon 
in the nucleus, the enhancement relative to the respective 
isoscalar operator tends to be slightly larger. In particu¬ 
lar, the constructive interference eigenvector for O^jO^ 
and OaIOq interference corresponds to WIMP-neutron 
spin-dependent scattering and is approximately a factor 
of 2 larger than the isoscalar O 4 rate. 


Since germanium, silicon, and xenon have similar prop¬ 
erties, the event rate in SuperCDMS and LZ is suppressed 
equally for most interference cases. However, there are 
a few notable exceptions. From Fig. 8 , we see that for a 
3GeV/c^ WIMP interacting via a pure isoscalar opera¬ 
tor, the event rate in SuperCDMS Si high-voltage tends 
to be at least an order of magnitude smaller than the 
rate in SuperCDMS Ge. When interference is consid¬ 
ered, the rate in silicon may become equal to or larger 
than that in germanium. As an example, the O 1 /O 3 
right-most destructive interference case in Fig. 9 corre¬ 
sponds to maximal Oi isospin violation in germanium 
{fn/fp —0.8) as discussed in [13]. For this choice of co¬ 
efficients, the rate in xenon and germanium is suppressed 
relative to pure isoscalar Oi scattering in that target by 
a factor of ^500 and ^ 2000 , respectively, while the rate 
in silicon is suppressed by a factor of ~100. A second 
instance of this suppression is seen for O^jOe interfer¬ 
ence at 3GeV/c^ in the second plot from the left: the 
rate in both silicon and germanium is suppressed, but 
the suppression in germanium is much larger, leading to 
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FIG. 9. Event rate suppression relative to 0\ scattering in LZ (black), SuperCDMS SNOLAB Ge iZIP (blue), and SuperCDMS 
SNOLAB Si (red) for interference in germanium, with interference ranging from constructive (left) to maximally destructive 
(right), as determined by the magnitude of the corresponding eigenvalue. The rate for SuperCDMS SNOLAB Ge high-voltage 
(light blue) is shown for the 3GeV/c^ case where LZ sees no events above threshold. The seven operator-operator interference 
cases are shown, as well as pure isoscalar 0 \, which is used as a reference point. 


a greater number of events observed in silicon. 

In addition, there exist several cases for higher WIMP 
masses where the rate in LZ is smaller than that in Super¬ 
CDMS Ge, despite LZ’s 100x larger exposure. Maximal 
destructive interference (right-most plot) for O^jO^ and 
0^1 Oq suppresses the event rate in xenon enough that 
SuperCDMS will see orders of magnitude more events 
than LZ, even for larger WIMP masses where LZ typi¬ 
cally has an advantage. For additional interference cases 
the rate in LZ is less than an order of magnitude larger 
than that in SuperCDMS Ge. Although the cases pre¬ 
sented here are arguably fine-tuned, the existence of re¬ 
gions of parameter space where interference suppresses 
the rate in one experiment by orders of magnitude rel¬ 
ative to another further supports the need for multiple 


experiments which use a variety of target elements. 


V. CONCLUSIONS 

The interaction between dark matter particles and nu¬ 
clei might be much more complicated than direct detec¬ 
tion experiments have typically assumed. The inclusion 
of new operators within the framework of an EFT might 
have profound consequences for current and proposed ex¬ 
periments. As a result, in this richer parameter space, 
data from multiple experiments with different targets is 
essential in order to determine the precise nature of the 
interaction. In addition, when modeling dark matter sig¬ 
nals, experiments must consider how an interaction due 
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to an arbitrary EFT operator can affect the energy dis¬ 
tribution of dark matter events. 

The importance of using multiple target elements to 
constrain dark matter interactions can already be seen 
when plotting limits from current experiments. As we 
have shown, the differences in target element properties 
lead to variations in the shape of the interaction strength 
versus mass limit curve. In addition, a combination of 
target elements can produce better constraints on dark 
matter, especially when considering multiple dark mat¬ 
ter interactions and the possibility of interference. This 
complementarity of different target elements will become 
increasingly important in the case of a statistically sig¬ 
nificant detection. 

The additional interactions introduced by the EFT for¬ 
malism become especially significant when experiments 
use statistical techniques which rely on assumptions 
about the shape of the dark matter recoil spectrum to dis¬ 
tinguish between background and a potential dark matter 
signal. Machine learning techniques, such as the boosted 
decision tree used in the SuperCDMS Soudan result [10], 
and likelihood analyses, such as the one performed on 
CDMS II low-energy data [34], require accurate mod¬ 
els of both the signal and the expected background. So 
far, direct detection experiments have focused primar¬ 
ily on building accurate models of their expected back¬ 
grounds, while assuming a simple signal model. However, 
mis-modeling the signal could also have significant con¬ 
sequences. If a WIMP signal that does not conform to 
the standard spin-independent assumptions is present in 
the data, it could produce unknown effects on the final 
result because it may not match either the signal or the 
background model. In the case of algorithms such as 
the optimum interval method that compare the observed 
events to the expected WIMP spectrum but do not at¬ 
tempt to subtract background, WIMP signal events may 


be interpreted as background, leading to limits that are 
too strict. 

These considerations become especially important as 
the community moves forward with the proposed G2 ex¬ 
periments. SuperCDMS SNOLAB and LZ will have un¬ 
precedented sensitivity to dark matter scattering for a 
wide range of WIMP masses, and the combination of 
target elements allows one experiment to verify a poten¬ 
tial signal seen by the other. However, the variation in 
signal strengths across EFT operators and experimental 
target elements could lead to experimental results that 
appear to be in conflict under the standard dark mat¬ 
ter assumptions. In particular, interference between op¬ 
erators can suppress the relative event rates by several 
orders of magnitude for germanium, silicon, and xenon. 
If the true dark matter interaction includes such inter¬ 
ference, it is possible that one experiment will observe 
a statistically significant signal while the other does not. 
The effective field theory framework can account for such 
apparent inconsistencies, and, in the event of a statisti¬ 
cally significant signal, it will pave the way for future 
likelihood analyses to determine the nature of the dark 
matter interactions. 
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Figure 1: Constructive interference eigenvectors for 4D O^jOg interference. (Published Fig. 6) 

In analyses undertaken after publication of the original article, a minor coding error (bug) has been identified in 
the utilized code. This bug is in the MATLAB code used to calculate detector target form factors and was identified 
by disagreement between the MATLAB code and the original Mathematica code [1]. The effects of this bug only 
occur in the interference between relevant operator pairs, such as operators 8 and 9 (Og/g). Any calculations using 
non-interfering pairs of operators or single operators are unaffected by the bug. The affected figures published in the 
original paper are reproduced in their corrected form here. The changes in these figures are relatively minor. 

The published Fig. 6 (Fig. 1 here) is one of the three affected plots. The directions of the constructive interference 
eigenvectors for each target in the 4D interference between operators 8 and 9 change, but the conclusion drawn from 
the figure does not: Fluorine is still the only target that is sensitive to the isoscalar and isovector components of 
operator 9. 

The published Fig. 7 (Fig. 2 here) is also affected by the coding bug. Because the direction of each constructive 
interference eigenvector changed, the point in the 4D parameter space at which this figure is made changed as well. 
For the new directions of the eigenvectors, the shapes of the corresponding event rates are very similar to those 
originally published with only minor changes. The main message of these plots does not change: For a given target’s 
constructive interference eigenvector, the other targets have lower or suppressed event rates in comparison to the 
chosen target. 

Finally, the published Fig. 9 is also affected. Due to its complexity and the CPU time required to reproduce it, this 
figure was not regenerated after correction. However, as with the other two figures, the changes in Fig. 9 would be 
minor and the main message remains unaffected: The event rates in different targets could be different by orders of 
magnitude depending on the possible constructive or destructive interferences among different operators. 


[1] N. Anand, A. L. Fitzpatrick, and W. C. Haxton, Phys. Rev. C 89, 065501 (2014), arXiv:1308.6288 [hep-phj. 
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Figure 2: Differential event rate evaluated at the 0%I0q constructive interference vector from Fig. 1 for fluorine 

(left) and germanium (right). (Published Fig. 7) 













